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In a previous paper, a binary black hole four-metric was presented in a post-Newtonian coro- 
tating coordinate system valid only up to the holes' apparent horizons. In this paper, I define an 
ingoing coordinate transformation that extends this corotating coordinate system through the holes' 
horizons and into their interiors. The motivation for using ingoing coordinates is that numerical 
simulations of black holes require the computational grid to extend inside the horizons. The coor- 
dinate transformation presented here makes the binary black hole four-metric suitable as a source 
of initial data for numerical simulations. 



I. INTRODUCTION 

In a previous paper an approximate solution to Einstein's equations representing two widely separated non- 
rotating black holes in a circular orbit was constructed by matching a post-Newtonian metric to two perturbed 
Schwarzschild metrics. The spacetime metric was presented in a single coordinate system valid up to the holes' appar- 
ent horizons. In this paper, I write the binary black hole four-metric from 1\ in coordinates that are corotating and 
post-Newtonian in the radiation and near zones and smoothly become ingoing near the black holes. This coordinate 
system is valid through the holes' horizons and covers the holes' interiors as well as the near and radiation zones. 
The metric components in this coordinate system are explicitly nonsingular on the black hole horizons. The metric 
presented here is promising as a source of initial data for numerical simulations of binary black holes. Since these 
simulations require the computational grid to extend inside the holes' horizons, the coordinate system used near the 
black holes in [l| is not suitable for numerical relativity. 

Let us begin with the metric near the first black hole, that is, in region I in the terminology of (see Fig. 1 
in 0). This metric is the Schwarzschild metric plus electric-type and magnetic-type tidal perturbations due to the 
second black hole, and is given in isotropic coordinates in Eq. (3.22) of The second black hole's tidal field rotates 
with angular velocity as seen by inertial observers in the first black hole's local asymptotic rest frame. However, 
the tidal perturbation's angular dependence — fit as given in Eq. (3.22) of 1] (I have replaced T in that equation 
with t) is singular at the first black hole's horizon. The reason is that the Schwarzschild time coordinate i, which 
is suitable for applying the technique of matched asymptotic expansions in the buffer zone around the black hole 
(see Q for details), is badly behaved at the horizon. Since our goal in this paper is to come up with coordinates 
valid through the horizon and inside the black hole, we must use a time coordinate T with the property that the 
hypersurfaces of constant time coincide with Schwarzschild time slices in the buffer zone but smoothly transition into 
ingoing Eddington-Finkelstein time slices which penetrate the horizon. The singular angular dependence — fit can 
be simply replaced by the nonsingular — fiT, with T as described above; this is discussed in further detail below. 

It turns out that isotropic coordinates are not a good starting point for an ingoing transformation. The analog of the 
ingoing Eddington-Finkelstein transformation, which is based on ingoing null geodesies of the Schwarzschild spacetime, 
is unsuccessful when applied to isotropic coordinates: the coordinate system remains singular at the horizon. Indeed, 
the isotropic radial coordinate is only defined outside the black hole. However, isotropic coordinates were used in 
to match a tidally perturbed Schwarzschild metric to the post-Newtonian near zone metric. It is therefore necessary 
to define a new radial coordinate that is equal to the (tidally distorted) isotropic radial coordinate in the buffer zone 
but transitions smoothly into the (tidally distorted) Schwarzschild radial coordinate near the black hole. 



II. INGOING COORDINATES 

Following the notation in [J, I denote the black holes' masses by mi and m2, and their coordinate separation in 
post-Newtonian harmonic coordinates by h. Let m = mi -I- m2, e = (m/&)i/2, and 17 = (1 — mim2 

/TO6)(m/63)i/2_ By 

assumption, e <C 1. 

Let us begin with the region I metric given in isotropic coordinates in Eq. (3.22) of 11|. Note that, in this paper, T 
and R denote the nonsingular time and radial coordinates described in Sec. ^ while in they denoted the isotropic 
time and radial coordinates. Set 17 = in Eq. (3.22) of Q and transform to Schwarzschild coordinates {i,f,9,(j)). 
This yields the metric g = gs -I- h; the Schwarzschild metric gg and the stationary tidal perturbation h are given in 
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Schwarzschild coordinates by 
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In this notation, dt, df, dO, and d(f> are coordinate one-forms and dt denotes the tensor product dt (g) dt. 

Let and (^2 be two numbers satisfying 2 < ^1 < ^2 < (6/toi)^/^. Define the new ingoing coordinates {T,R,6,( 
by 



t = T — 2mi In 
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f = i? + mi(l + ^)r;(i?). 
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The functions ^/'(-R) and r/(i?) must satisfy the following properties: (i) for R > C2W1, ip{R) — and ri{R) = 1 
so the coordinates are isotropic; (ii) for R < Ci™!, ^{R) = 1 and ri{R) = so the time coordinate is Eddington- 
Finkelstcin and the radial coordinate is Schwarzschild; and (iii) for Ci"^i < R < C2mi, tp{R) and ri{R) smoothly 
and monotonically vary between their constant values outside this interval. The transition points R = ^imi and 
R — C2fni and the functions ip{R) and ri(R) can be chosen freely as long as the above properties are satisfied. Since 
df/dR = 1 + (mi + rnl/m)ri'{R) - mlr]{R) / iR"^ > 1 - ml/iR'^, we have df/dR > for > 2mi, and so Eqs. ® 
and I0J define a valid coordinate transformation. 

Transforming and j^Jl using Q and I0J, we obtain 
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and f is given in terms of R in Since the linearized Einstein equation is a tensor equation (see, e.g., Eq. (7.5.15) in 
0) and we have only performed a coordinate transformation, h remains a solution to this equation on the background 
gs. However, the 0-dependence of the perturbation h does not correspond to the second black hole's rotating tidal 
field — we set = above. To remedy this, we simply change to — OT in the components of h, which yields a 
new perturbation h. Note that this replacement is not a coordinate transformation; a new tensor h is defined. Also 
note that T becomes an ingoing coordinate near the horizon, so the time dependence (j) ~ causes no problems 
on the horizon. This simple remedy works for the following reason: if we solve the linearized Einstein Eq. order by 
order in e, then time derivatives of the components of h produce factors of niiil ^ c' and can thus be neglected. The 
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perturbation h is given by 
h = 
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We have now arrived at a metric g = g5 + h, where gs is given in Q and h in (|10|l . which is vaHd from the black 
hole's interior up into and through the buffer zone around the hole, and is written in coordinates that are well behaved 
throughout this region. 

III. TRANSFORMATION TO COROTATING COORDINATES 

The next step is to transform the metric gto corotating post-Newtonian coordinates {t,x,y,z) using the trans- 
formation given in Eqs. (4.22) and (4.23) of Q. This transformation contains a rotation that can be performed by 
first defining (p = cj) — flT and then setting F = i?sin0cos(/?, A = i? sin sin 1^9, and Z = RcosO. To complete the 
transformation, define the functions Pap{x,y, z) for a,f3 = 0, ..,3 to be components of the metric g in coordinates 
(T, r,A, Z); write these components as functions of {x,y,z) using Eqs. (4.22) and (4.23) in T|. The functions Pa/s 
are given by 
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where e = {m/by/^, Q ~ {1 — mim2/mb){m/b^y/^ , R — (T'^ + A2 + ^2)^/2, and F, A, and Z are to be expressed in 
terms of {x,y,z) via Eqs. (4.22) and (4.23) in 1]. The functions f{R), ^'(i?), H{R), and f{R) are given in Eqs. 
O, ©, and ®. 

The remainder of the coordinate transformation from black hole coordinates (T, R, 9, (/)) to corotating post- 
Newtonian coordinates (t, y, z) can be done exactly as in [J. The final metric in region I (see Fig. 1 in 1]) is 
given by Eq. (4.27) of but with Pap taken from Eq. <fTT|l above. Define Paf3 to be Paf3 (as given in l lTT|l ') with 
TOi and m,2 exchanged. Then the metric in region II (that is, near the second black hole; see Fig. 1 in [l|) is given 
by Eq. (4.28) in 1], but with P^^ taken from here. Note that the final metric components are everywhere explicitly 
independent of time t. 

To summarize, the expressions for the metric components given in Sec. V of '1] are valid in ingoing coordinates 
if Pap (and Pap) are taken from Eq. I^ll]) and the functions r(R), "^'{R), H{R), and f{R) are taken from Eqs. 
Q), 0), and with ip{R) and r](R) having properties (i)-(iii) given below Eq. In order to show explicitly that 
these components are nonsingular at the black holes' horizons, I write out P^p for R < Ci^i- In this region, ip = 1, 
^{R) = 2mii?-i(l - 2mi/R)-\ 77 = 0, = 1, f = i?, and f{R) = 1 - 2mi/R. Therefore, for R < Cimi, F„/3 are 
given by 
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P23 = P32 = + -^z{2z' - R') + -^(3r' - (12) 

where e = (m/fe)^/^, = (1 - mim2/mb){m/b^)^/'^ , R — (F^ + + Z^)^/^, and F, A, and Z are to be expressed 
in terms of {x,y,z) via Eqs. (4.22) and (4.23) in H. Note that the quantities Pap in Eq. (|12|l are all finite at the 
horizon R = 2mi. 
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